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By Harry L. Runyan, Donald S. Woolston, 
and A . Gerald Rainey 


SUMMARY 


A recently published analytical investigation of the effects of 
wind-tunnel walls on air forces on an oscillating wing at subsonic speeds 
(NACA TN 2552) demonstrated the possibility, under certain conditions, of 
the existence of large tunnel-wall effects, associated with an acoustic 
resonance phenomenon. In the present paper the integral equation defining 
the problem is treated further and is presented in a form adapted to 
calculations. Application is made to the case of a pitching airfoil in 
a given wind tunnel at a Mach number of 0.7 for frequencies ranging from 
zero to beyond the frequency of acoustic resonance. The calculated lift 
and phase angle are compared with the results of experimental measurements 
made in the Langley 4.5-foot flutter research tunnel. Good agreement is 
obtained between the calculated and the experimentally determined magni- 
tr of the lift and phase angles. 


INTRODUCTION 


In be evaluation of results obtained by measurement of the air 
forces on an oscillating wing in a wind tunnel, the question of the 
effect of the tunnel walls arises. A theoretical treatment of the effect 
of wind-tunnel walls for the case of two-dimensional incompressible flow • 
has been made by several investigators and reported in references 1, 2, 
and 3- In these papers the influence of the tunnel walls is found to be 
comparatively small for most cases. Extension of the treatment of the 
problem to include the effects of compressibility of .the fluid has been 
reported in reference 4. It is shown in this reference that, for certain 
conditions of wing frequency, tunnel height, and Mach number, the oscil- 
lating wing could excite the' air in the tunnel to the point of resonance 
and that, in the absence of damping, infinite transverse velocities could 
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exist. In the neighborhood of such a resonant condition one would 
expect very large wall effects and it would be extremely difficult to 
approximate the free-air condition. 

In reference 4 sin integral equation is derived which relates the 
downwash distribution to the lift distribution for an airfoil of infinite 
aspect ratio oscillating in a wind tunnel. A main purpose of the present 
paper is to present a method of solution of this equation and to demon- 
strate some applications. The method of solution is one of collocation, 
similar to that used by Possio (ref. 5) and Frazer (ref. 6) , in which 
the downwash is satisfied at selected control points along the wing 
chord. Accuracy of the solution is increased as the number of control 
points is increased, and in the present investigation successive calcu- 
lations based on one, two, and three control points have been made - 

Another objective of this paper is to make comparisons with experi- 
mental results and for this purpose measurements have been made of the 
magnitude and phase relationships of the lift on a two-dimensional wing 
oscillating in pitch. Measurements were made at a constant Mach number 
(0.7) for a range of frequencies of oscillation. These results are com- 
pared with the results of calculations made for the same conditions. 

As stated in reference 4, the integral equation for the case of no 
tunnel walls checks the results of Possio. For the case with t unn el 
walls and for the limiting steady-flow case of zero frequency, it is 
possible to obtain a mathematical check with some existing results; this 
is shown in an appendix. 


SYMBOLS 


A n coefficients in series expression for lift distribution 

(eq. 11) 

a axis of rotation measured from midchord, positive rearward, 

based on half-chord 

b wing half-chord, ft 

c velocity of sound, ft/sec 

h displacement of wing in vertical translation 

H height of tunnel 

( 2 ^ 

H n v } Hankel function of the second kind 
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k reduced frequency, boo/u 

K(M,z) + K(M,z,H) kernel of integral equation 
L(x 0 ), L(0 O ) lift distribution, lb/ft/unit span 
M Mach number, U/c 



U stream velocity in chordwise direction, ft/sec 

w(x) vertical induced velocity (perpendicular to chord), ft/sec 

,x, x Q , | chordwise coordinates 
y vertical coordinate 

z = k( x — x Q ) 

a angular displacement of wing in pitch 



to circular frequency of oscillation, radians/sec 

u> res circular frequency at resonance, radians/sec 

p fluid density, slugs/ft3'. 

Subscripts: • 

n 0, 1, 2, . . . 

Following equation (6) variables i, x, x 0 , H, and y are con- 
sidered as nondimens ional length quantities based on the half-chord b. 
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This section is given in two parts: First the integral equation is given in a form suit- 

able for computation, and second the method of calculation is shown. The basic assumptions in 
the derivation of the integral equation of reference 4 were the usual ones of linearization, 
nonviscous fluid, and thin bodies. In addition, it was assumed that the tunnel walls form 
perfect reflecting surfaces for the disturbances. 
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The dimensional form of the equations is retained, but later a nondimens ional form is given. 
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in which use has been made of 



Although the kernels K(M,z) and K(M,z,H) are similar, it may 
be noted that K(M,z) is singular since, at x = x Q , R 0 is zero,- and 

for small argument, H 0 ^(x) behaves as log e x and behaves 

as l/x. The argument R n of the Hankel functions in K(M,z,H) does 
not approach zero and therefore the terms in- the bracket for this part 
of the kernel are not singular. It will be seen later, however, that 
the series representation of K(M,z,H) becomes large as resonant fre- 
quencies are approached and is infinite at the resonant condition. 

The integral equation for the case of no tunnel walls checks the 
results of Possio. For the case with tunnel walls and for the limiting 
steady-flow case of zero frequency, it is possible to obtain a mathe- 
matical check with some existing results; this is shown in an appendix. 


Method of Solution 

A method of using equation (7) to determine the aerodynamic forces 
on a wing oscillating in the presence of tunnel walls is now discussed. 
The method under consideration is one of collocation, similar to that 
used in references 5 and 6 for the case of no tunnel walls. The method 
involves writing the expression for the downwash at selected control 
points, that is, reducing equation (7) to a set of linear algebraic 
equations. The unknowns in the set of simultaneous equations are the 
various coefficients of the expression which is assumed for the loading. 

This expression which is assumed for the present case is a trigono- 
metric series expansion which satisfies the Kutta condition at the 
trailing edge and which has the proper type of infinity at the leading 
edge. This expression is 

L ( x \ 6 oo 

_L_2i = Ao cot - 2 - + A n sin n0 o = L(0 O ) (10) 

P u n=l 
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where x Q = -cos 0 O and the A n 's are unknown coefficients to he deter- 
mined in accordance with the downwash w(x), which is known from the 
motion of the wing. It is desirable to rewrite equation ( 7 ) in terms of 
the new variable Q 0 , and equation (7) then assumes the form 



r\ ^ 

f ^ 

w(x) = Uk 

/ L(0 o )K(M,z) sin 9 0 d6 0 + 

/ L(0 o )K(M, z,H) sin 9 0 &9 0 


Jo 

Jo 


= Ukt 1 + I 11 ) (11) 

The integral I 1 of equation (ll) to be discussed, first is the integral 
expression as first derived by Possio for the condition of no tunnel 
walls and has been treated by several investigators. As pointed out 

previously, the kernel K(M,z) of I 1 is singular at x = x 0 . Since, 

in handling the complete integral I*, numerical or graphical means are 
necessary, it is convenient to isolate the singularity and express it 
in an integrable form. This has been done by Schwarz (ref. 7) who 
expanded K(M, z) in the vicinity of z = 0 in the form 


K(M,z) 


F(M) + 
z 


iG(M)log e 


z l 


+ K x (M,z) 


where F(m) = G(M) = — , and K-i(M,z) is a tabulated function 

which is no longer singular. Substitution of this expression in the 
integral I 1 of equation (ll) gives 


I 1 = / L(0 O ) sin 0 o K 1 (M,z)d0 o + 

Jo 


L(0 O ) sin e 0 


F(M) 


■H 


+ iG(M)log e |z| d0 o (12) 
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With the use of equation (10), the integral of equation (12) con- 
taining K]_(M,z) becomes 



L(0 o )sin 9 q KpM,z)d0 o 



(1 + cos 0 o )KpM,z)d0 o + 


oo p Tt i 

^ A n / sin 0 O sin n0 o K-^M, z) d0 o ( 13 ) 
n=l O 0 


By utilizing the values of KpM, z) tabulated by Schwarz, the indicated 
integrations of equation (13) may be readily carried out by numerical or 
graphical means. 

Putting z = -k(cos 0 - cos 0 O ) gives for the second integral 
of in equation (12) 


L(0 o )sin 0 O 


+ iG(M)log e | z | 




1 . L(0 o )sin 0 O 

iF(M) / — — d0 n + 

k Jq cos 0 - cos 0 O ° 


iG(M) J L(0 o )sin 0 O log e | cos 0 - cos 0 o |d0 o + 


iG(M)log e k / L(0 o )sin 0 O d 0 O 

J 0 


= -i F(M)I 0 + iG(M) Ip + iG(M)l 2 log k 


. ( 14 ) 


The integrals I Q , Ip, and I 2 have been evaluated by Frazer (ref. 6) 
and have the following values: 

00 

I Q = -JtA 0 + jt ^ A n cos n0 (15a) 

n=l 
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II = -it^A 0 + i A^log e 2 - nA 0 cos 0 + Ap cos 26 


lV A 

p Z H n 

n=2 


cos(n + l)0 cos(n - 1)0 


n + 1 


n - 1 


( 1513 ) 


I 2 - *(a 0 + i A^ 


(15c) 


With the use of these expressions in equation (l4) and with the results 

of the integration of equation (13), the- integral I 1 of equation ( 11 ) 
is now expressed solely in terms of the unknown coefficients A n . It 

remains now to obtain I • of equation (ll) in a similar form. 

The kernel K(M,z,H) of I 1 " 1 ", given in equation (9), is the sum 
of four infinite series; namely 


K(M,z,H) = 


,-iz 


2(3 


1Z 


-e 


p XI (-D n H o ( 2 )0iR n ) + 

n=l 


l» . J 


a (-D" + 

n=l \ P / 


-ik| 


ik X(-D n / e p2 H o (2) (p^ 2 +.P 2 (nH) 2 jd£ + 


n=l u 0 


iz 



n=l 


(_l) n iM(x - x Q ) 
R n 


Hi (2) (pR n ) 


p -iz 

= %3~( c l s l + C 2 S 2 


* C 3 S 3 


+ VO 


( 16 ) 


where the S n 's denote the indicated infinite summations and the C n r s 
the respective multipliers. 
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The series Sp and Sp of equation (l 6 ), which involve the Hankel 

(p) 

function H 0 V , may be put in a more rapidly convergent form according 
to reference 8 by introducing the variables p and e, where 

p = m. 

P 2itp 


and 


x - x r 


e = 


3H 


Accordingly, the series Sp and S 2 can be written as 

CO 

Si' -51 (-D n H 0 ( 2 )(pR n ) 

n=l 


n=l 


(-l) n H 0 ( 2 ) ( 2 np \/e 2 + n 2 ) 


exp(-it€^T- 4p 2 ) _ Tf- . 

exp 

_-ne y(2n + l) 2 - Up 2 

VlT l,p2 n=l 

\/(2n + l) 2 - 4p 2 


exp 




2n - l) 2 - 4p^_. 


\f( 2 n - l ) 2 - %> 2 


l- jrH D ( 2 )(2nep) 


(17) 


and 


s 2 - 51 (-!)■> 

n=l V P / 

H 0 ( 2 )(2nnp) 


n=l 


-n + 2i (7 + log e 2p) + 4i Y~ ( - 1 i \ 

L - !)2 - 4 p 2 211 - 7 


(18) 


where Eulers' constant 7 = 0.577215. 
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Series S^ was evaluated in a direct manner by employing tables 

of the Hankel function and by using for large values of the argument 
the approximation 



With the aid of these series S]_, S 2 , S 3 , and S^, the kernel 

K(M,z,H) may be evaluated. The integral I 11 of equation (ll) maybe 
expressed in terms of the unknown coefficients A n upon substitution 
of the assumed expression for L(0 O ) and by graphical integration. 

The integrals of equation (ll) are determined in this manner for a 
selected number of control points and equated to the expression for the 
downwash . 

The expression relating the downwash to the motion of a plunging 
and pitching wing is 


w(x) = h + Ua + b(x - a)a ( 21 ) 

With the assumption of harmonic motion, equation (2l) is used to calcu- 
late w(x) for values of x appropriate to each of the control points. 
A set of simultaneous equations can then be written, the number of which 
corresponds to the number of control points assumed and to the number 
of terms retained in the series for L(0 O ) . The unknown coefficients 
may now be determined by solving these simultaneous equations. It may 
then be shown (see, for instance, ref. 9 ) that the total lift and moment 
are given by 


-L 

jtpbU 2 


= i{ 


v + 



> 


Mq 

:tpb 2 U 2 




( 22 ) 


where is the moment about the midchord. 
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Summary of procedure for calculating tunnel-wall effects .- As an 
aid in calculating the effects of wind-tunnel walls on an oscillating 
wing the following summary of the procedure is given. It is assumed 
that the Mach number, speed of sound, wing chord, frequency of oscil- 
lation, and tunnel height are known. 

(1) Select the position x of each of the control points for which 

equation ( ll) will be written. First, treat I* of equation (ll) as 
given by equation (12). 

( 2 ) For each value of x, evaluate the first integral on the right 
side of equation (12), in the form given by equation ( 13 ) , by numerical 
or graphical means. 

( 3 ) Evaluate the second integral of equation ( 12 ) in the form given 
by equation (lb) by use of equation (15). Add the results to the results 

of step (2) to obtain 1^ in terms of the unknowns A n . 

( 4 ) For I 11 of equation (ll), first evaluate K(M, z,H)- given by 
equation ( l6) by use of the series given in equations (17) to (20). 

( 5 ) For each value of x substitute the values of K(M,z,H) from 
step ( 4 ) in I 11 of equation (ll) together with L( 0 O ) from equa- 
tion (10) and perform graphical or numerical integration to give 

in terms of the A n . Combine these results with the results of step (3) 
for the corresponding values of x. 

(6) From equation ( 2 l), with the assumption of harmonic motion, 

calculate w(x) for values of x appropriate to each of the control 
points. ' 1 

( 7 ) Equate results of steps (5) and (6) for corresponding values 
of x to obtain a set of simultaneous equations in the unknowns A n . 

(8) Solve the equations from step (. 7 ) for the A n 's. 

( 9 ) Determine lift and moment from equations ( 22 ). 

Application to a specific problem .- The foregoing theory has been 
applied to the particular case of a wing oscillating in pitch about the 
midchord (a = 0 , h=0 in eq. (21)). Calculations have been performed 

for M = 0 - 7 > H = 3 *802 feet, b = i foot and c = 531 feet per second 

for various frequencies of oscillation from 0 to 60 cycles per second. 
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In performing the calculations, the number of, terms of the series 
for lift distribution (eq. (10)) and thus the number of control points 
required to obtain satisfactory accuracy was not known. Accordingly, 
the calculations were performed by increasing the number of terms of 
the series until the solutions were in reasonable agreement. First, 

the initial term of the series for lift distribution A 0 cot ^ was 


used together with a control point at ^ - chord position. Then two terms 

e Q 4 

A 0 cot — + Ai sin 9 0 and two control points located at 1- and - -chord 


positions were used. Finally, the first three terms of the series and 
control points at and ^ - chord positions were utilized. In 

calculations were made for the case of one doublet placed at 

the i- chord position and for the downwash satisfied at the ^ - chord 

. 4 

position. The results of these calculations are shown in figure 1. 


In the calculations of the series for the determination of K(M,z,H) 
it was found that as the resonant condition was approached, more terms 
were required for convergence. In all cases, however, a sufficient num- 
ber of terms of the series was taken to assure convergence in the third 
decimal place. Series was the most slowly convergent and in one 

case required 240 terms in order to assure convergence. 


EXPERIMENT 


The experimental part of the investigation of the effect of wind- 
tunnel walls was conducted in the Langley 4 *5— foot flutter research 
tunnel in which a 2-foot by 4-foot rectangular section has been tempo- 
rarily installed for testing two-dimensional models. This tunnel is of 
the closed- throat, single-return type employing either air or Freon-12 
as a testing medium at pressures from 1 atmosphere down to about 
1/8 atmosphere. Since Freon-12 has a speed of sound equal to about 
one-half that of air and since the critical tunnel frequency varies 
directly as the speed of sound, the experiments to be discussed were 
conducted in Freon-12 so that the resonant frequency could be surveyed. 
The Mach number of the tests was M = 0.7 and the Reynolds number 

was 4.01 X 10 6 ., The frequency range was from 0 to 57-54 cycles per 
second, and the amplitude was about 2.4°. 

The model was of 1-foot chord having NACA 65-010 airfoil sections 
and completely spanned the 2-foot dimension of the test section. The 
gaps between the wing and the tunnel wall were sealed by end plates 
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which rotated with the model. The wing was oscillated in pitch about 
the midchord by a direct-drive eccentric-cam system powered by an induc- 
tion motor with variable frequency supply. 

The lift on the wing was obtained by electrical integration (ref. 10) 
of the outputs of 12 model 49-TP NACA miniature electrical pressure gages 
(ref. 11) arranged to indicate the differential pressures between orifices 
on the upper and lower surfaces at the midspan position. The angular 
motion of the midspan position was indicated by resistance wire strain 
gages attached to a torque rod running through the center of the hollow 
wing. One end of the torque rod was fixed to the center of the wing and 
the other end was fixed to the wind-tunnel wall. 

The amplitude of the fundamental frequency of both the lift and 
position were indicated on an alternating- current vacuum-tube voltmeter 
attached to the output of a variable -frequency, narrow -band-pass filter. 
The phase angle between lift and position was determined by use of a 
special electronic timing system. 


RESULTS AND DISCUSSION 


A comparison of analytical and experimental results obtained for 
the lift on a wing pitching about its midchord is shown in figures 1 
and 2. Variations in the magnitude of the lift and in phase angle are 
given as functions of the frequency of the pitching oscillation for a 
constant value of Mach number. 

The ordinate in figure 1, the magnitude of the lift, is presented 
as a ratio of the value obtained either experimentally or theoretically, 
with consideration of the tunnel walls, to the theoretical value for the 
condition of no tunnel walls. As previously mentioned, in investigating 
the number of control points required to give sufficient accuracy, a 
series of calculations was made at each pitching . frequency with each 
succeeding calculation including an added control point. The result of 
each of these calculations is shown in figure 1. 

The abscissa in figure 1 is the ratio of the frequency of the 
pitching oscillation to a frequency calculated for the resonant condi- 
tion. Since in the experiment the velocity of sound of the mediums in 
which the tests were made differed slightly from that considered in the 
theoretical work, the calculated resonant frequencies were, not quite the 
same. The resonant frequency for conditions of the calculations was 
49.5 cycles per second; whereas for conditions of the experiment, the 
resonant frequency was 52.6 cycles per second. The pitching frequencies 
for the calculated and the experimental results were divided by the 
appropriate value of the resonant frequency. 
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The experimental results show the existence of large tunnel-wall 
effects and demonstrate the existence of a resonant phenomenon as pre- 
dicted in reference 4 . The theoretical results are in very good agree- 
ment with experiment with regard to predicting the value of the resonant 
frequency. At low values of the frequency ratio the theoretical results 
for the lift ratio given in figure 1 indicate good agreement between the 
calculations for 1, 2, and 3 control points and are in good agreement 
with the experiment. At values of the frequency ratio nearer unity there 
is a wider divergence between the various solutions and between theory 
and experiment. This indicates that near res.onance additional control 
points or possibly different control-point locations would be required 
to give a more accurate approach to the resonant condition. 

As a matter of interest, and as a partial check on the mechanics 
of the procedure, calculations of the aerodynamic lift and moment coef- 
ficients for. the condition of no tunnel walls were performed with three 
control points, and excellent agreement was obtained with the tabulated 
results of Dietze. 

In figure 2 the phase angle between lift force and position of the 
wing is plotted against the same frequency ratio <^/^ res used in fig- 
ure 1. Although the variation of phase angle with frequency in free air 
is almost linear, the phase angles found both experimentally and by the 
theory including wall effects show a very large drop as the resonant 
frequency is approached and change abruptly. The calculated values of 
the phase angle are in good quantitative agreement with those measured 
experimentally. 

' Since the underlying theory of an oscillating wing in a compressible 
fluid is one of acoustics, it seems reasonable to conjecture that some 
application of the classic principles of sound absorption to the test 
section of a tunnel might be a method of eliminating or greatly reducing 
this new type of tunnel-wall interference. Any method of sound absorp- 
tion employed, of course, should not interfere with the basic flow in 
the tunnel. One tunnel configuration, which has been suggested for the 
steady-flow case and which may be desirable for the present case, is the 
ventilated tunnel. 


CONCLUSIONS 


This paper presents the results of an application of the theory as 
developed in NACA TN 2552 (which treats the' effect of the presence of 
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tunnel walls for oscillating two-dimensional compressible flow) and com- 
pares the results with some experiments at a Mach number of 0.7. The 
following conclusions may he enumerated: 

1. Some experimental results of lift and phase, angles made at a 
constant Mach number and for various oscillating frequencies show the 
existence of large tunnel-wall effects and demonstrate the existence of 
a resonant phenomenon as predicted in NACA TN 2552. 

2. A method of calculating the effect of tunnel walls on the 
oscillating two-dimensional air forces is presented. 

3. Comparison of the experimental lift with the theoretical lift 
which includes the effect of tunnel walls indicates good agreement 
throughout the frequency range. 

4. A comparison of the phase angles as determined by experiment and 
as calculated by the theory which includes the effect of tunnel walls 
shows good agreement throughout most of the frequency range and shows a 
very large phase-angle change as the resonant frequency is approached. 

It is observed, however, that the theory without tunnel walls does not 
give either the trend or the magnitude of the experimental results. 

5* At the Mach number of these tests (0-7) the effects of the tun-: 
nel walls is very great at the resonant frequency and the large effect 
extends over a considerable range of frequencies. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics , 
Iangley Field, Va. 
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APPENDIX 

REDUCTION OF INTEGRAL EQUATION TO THE CASE OF ZERO FREQUENCY 


In this appendix, the integral equation for the downwash for a wing 
oscillating in a compressible medium and enclosed by wind-tunnel walls 
is reduced to the zero-frequency condition. 

If equation (7) is written as 

w(x) = lim -i- f L(x 0 )HIk.(M,z) + ccK(M,z,H)ldx 0 (Al) 
oo ^ 0 pU^ J -*b — 1 


and the limit is taken as a> — > 0, it will be found that all the terms 
of K(M,z) and K(M,z,H) are canceled qxcept terms involving . 

These terms become infinite, however, as u) — ^ 0 so that the asymptotic 
expansion for very small values of the argument may be used; therefore. 


and 


H! (2) (|!Rn) ~ 


2 


1Z 


lim a f > 2 E^) - - 2Mce (* - *°> 

^ * 0 n * (X - X 0 ) 2 + p 2 (nH) 2 


The vertical induced velocity may then be written as 
>b 


w(x) = / L(x 0 ) 

2«pU 2 J _b 


^- + 2 jK-i) n 


X - X, 


x - x_ 


n=l (x - x 0 ) 2 + 3 2 (nH) 2 


dx. 


(A2) 


or 


w(x) = / L(x 0 ) 

2pU 2 H J ' 


00 n ( X - X o) 

+2^ (-l) n m 

m ( x ~ x o) n=l 


ffg ( x - x o) 2 . 2 2 

3% 2 


dx r 


(A3) 
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This equation may "be written as 



The additional induced velocity due to the presence of tunnel walls 
for the steady case in compressible flow is given by equation ( ho ) of 
reference 12. Equation (A2) can be reduced to the same form by making 
the approximation that the airfoil chord is small compared to the tunnel 
height. 
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Figure 1.- Plot of ratio of magnitude of experimental and theoretical 
lift, including tunnel-wall effects, to the theoretical lift without 
tunnel-wall effects against frequency ratio of oscillation to resonant 
frequency. M = 0.7. 








